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ABSTRACT 

Context. Dust emission at sub-millimeter to centimeter wavelengths is often simply the Rayleigh- Jeans tail of dust particles at thermal 
equilibrium and is used as a cold mass tracer in various environments including nearby galaxies. However, well-sampled spectral 
energy distributions of the nearby, star-forming Magellanic Clouds have a pronounced (sub-)millimeter excess (Israel et al., 2010). 
Aims. This study attempts to confirm the existence of such a millimeter excess above expected dust, free-free and synchrotron emission 
and to explore different possibilities for its origin. 

Methods. We model near-infrared to radio spectral energy distributions of the Magellanic Clouds with dust, free-free and synchrotron 
emission. A millimeter excess emission is confirmed above these components and its spectral shape and intensity are analysed in light 
of different scenarios: very cold dust, Cosmic Microwave Background (CMB) fluctuations, a change of the dust spectral index and 
spinning dust emission. 

Results. We show that very cold dust or CMB fluctuations are very unlikely explanations for the observed excess in these two 
galaxies. The excess in the Large Magellanic Cloud can be satisfactorily explained either by a change of the spectral index due to 
intrinsic properties of amorphous grains, or by spinning dust emission. In the Small Magellanic Cloud however, due to the importance 
of the excess, the dust grain model including TLS/DCD effects cannot reproduce the observed emission in a simple way. A possible 
solution was achieved with spinning dust emission, but many assumptions on the physical state of the interstellar medium had to be 
made. 

Conclusions. Further studies, using higher resolution data from Planck and Herschel, are needed to probe the origin of this observed 
submm-cm excess more definitely. Our study shows that the different possible origins will be best distinguished where the excess is 
the highest, as is the case in the Small Magellanic Cloud. 

Key words. Magellanic Clouds, SubmillimeterTSM, Radio continuum: ISM, ISM:dust, GalaxiesTSM 



r'"] 1. Introduction ified black-body (S v ^BvOTdusi) where f3 is frequently taken 

as ~ 2). Because of the wavelength dependence of thermal dust 

»H ; (Sub-)millimeter to centimeter emission in galaxies is thought emission? m iiii me ter dust emission is often used as a tracer of 

to be quite simple, originating from a combination of dust emis- ±e c old molecular ga s res ervoir in a galaxy (e.g.lGuelin et all 

sion, free-free and synchrotron emission. Thermal free-free ra- [T993I: iDunne & Eales 200ll: lRubio et alll2TO[WeufetliT2 008: 

diation originates in the ionized gas of Hn regions and its emis- |yi aha Ms et al.l2008t This information is particularly interesting 

sion is nearly flat in this optically thin regime of the spec- since cold Ui is almost im p OSsible to observe and the use of CO 

trum (S v oc v- 01 ). Synchrotron radiation is emitted by rela- i ine e mission is not straightforward, in particular at l ow metallic- 

tivistic electrons accelerated in magnetic fields and its spec- it ies (IIsraelll988tlLeaueux et al Jl994t llsraell 1997bh . Millimeter 

trum is also a power law (S v oc v "™« )■ The spectral index of dust emission cou i d be a good alternative to the CO molecule 

synchrotron emission and the relative contribution of these two f or tracing cold matter, also in low metallicity environments 

emission mechanisms to the intensity observed at a given wave- (Thronson 1988l lIsraefll9973lRubio et : al.l2004tiBot et .8112 007. 

length is constrained by radio observations. Dust emission ob- |2Q1Q|) 
served at far-infrared and millimeter wavelengths is produced 

by large dust grains in thermal equilibrium that emit as a mod- This simple picture is however thrown into doubt by obser- 
vations of the Milky Way and nearby galaxies. Several studies 
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trace large amounts of very cold dust with temperatures below 
10K (Reach etall 119951: IChini etal.lll995t iKriigel etail 11998" 



Pope scu et al.l l2002t iGalliano et all 120031: iDumke et alj 1200- 
But, as already noted by iReach et aH (1 19951) . the spatial cor- 
relation of this very cold dust with the warm component in 
our galaxy argues against this interpretation. In nearby galax- 
ies also, the scenario of ver y cold dust is rejecte d based on 
different physical arguments dLisenfeld et al.l 120021 e.g.). This 
sub-millimeter excess is then alternatively attributed to varia- 
tions of optical properties of dust grains. Fo r example, low spec- 
tral emissivity index values, /3, might apply dPupac et al.l l2003b; 

Variations of the spe ctral emissivity in- 



Lvity i 
lAguirre etall 120031) . 

dex with dust temperature have been noted (Dupac et al. 2003a; 



IDesert etaill2008l) . These variatio ns could be relate d to physi- 
cal properties of amorphous solids dMenv et al.l20 07). However, 
in nearby galaxies, such variations in dust properties might be 
difficult to disentangle from the effect of mixing different dust 
components, with different dust temperatures, along the line of 
sight. Data at numerous wavelengths covering the far-infrared 
(FIR) to the millimeter regime are required to better constrain 
the emission from dust in nearby galaxies. 

At longer wavelengths (~ 30GHz, ~ 1cm), another 
excess of emission has been ob served in our galaxy and 
is called "anomalous dust" (e.g . Draine & Lazarian 1998a; 
Lagache 2003; de Oliveira-Costa et al.l l2004t [F inkbeiner et al] 
2004 ICasassus et all 120061: iMiville-Deschenes et all 12008: 
Dobler & Finkbeinen 120081) . This excess is best explained by 
small rotating dust grains ("spinning dust" iDraine & Lazariarj 
1998b). Anomalous dust emission remains poorly constrained 
by observations (in particular its variations with the environ- 
ment). Theoretical models predict that the peak frequency of 
spinning dust emission should shift with grain size or density 
dAli-Haimoud et al.l l2009t lYsard & Verstraetd l2010l) . What is 
the influence of spinning dust emission on submillimeter to radio 
spectral energy distribution of galaxies? 

The Small and the Large Magellanic Clouds (SMC and 
LMC) are two of the nearest galaxies to ours (located at dis- 
tances of 60 and 50kpc respectiv e ly ([D eb & Singh 2010; 
ISzewczvketaLll2009t ISchaefedl2008l: ISzewczvk et alJ l2008)). 
As such, they have been extensively observed in all wavelength 
regimes. Their low metall icities (~ 1/ 2 and 1/6 for the LMC 
and the SMC respectively (Pagel 2003)) provide an opportunity 
to test our understanding of astrophysical processes in different 
physical conditions from our galaxy. 

Ilsrael et ail d2010t) build spectral energy distributions (SEDs) 
of the Magellanic Clouds by combining selected literature flux 
densities with COBE and WMAP data. The spatial resolution is 
low (~ 1" beam), but the resulting SEDs have an unprecedented 
sampling of wavelengths, particularly at the far-infrared to the 
radio regime and show a pronounced (sub-)millimeter excess. In 
this papier, we analyse these SEDs of the full SMC and LMC to 
assess the contribution of dust, free-free and synchrotron emis- 
sion, in order to quantify and study this excess emission and its 
possible origin. Our study is a prototype of future work on more 
distant galaxies that will use Herschel (60-600^m) and Planck ( 
up to ~ 10 mm) data. 

In section [2] we present a fit of the integrated spectral en- 
ergy distributions of the Magellanic Clouds, at infrared and ra- 
dio wavelengths, with dust, free-free and synchrotron emission. 
Above these three components, a significant emission excess 
is observed in the millimeter-centimeter domain. We explore 
whether this excess could be associated with very cold dust (Sec. 
@J, with cosmic microwave background fluctuations (Sec. |3), 



with the amorphous nature of dust grains (Sec. [5]>, or with spin- 
ning dust emission (Sec |6j. 



2. Integrated spectra of the Magellanic Clouds 

We used t he sp ectral energy distribution determined by 
Israe l et al. I (12010) for the Small and the Large Magellanic 
Clouds. These integrated spectral energy distributions are shown 
in Figure Q] and the in frared and radio em ission is fitted with 
three components: the IDraine & Lil (120071) dust model, a free- 
free component and synchrotron emission. 

The choice of a model to repr esent dust in the 



Magellanic C louds i s still problematic (e .g. Meixner et al. | 20 1 
Gor don et al.l 120101: iRoman-Duval et al.ll2010t iKim et aLlpOl 
Honv et al.ll2010h . We chose to use the lDraine & Lil d2007l) dust 
model, which is well suited to studies of galaxies spectral energy 
distribution and has been applied to a number of n earby galax- 
ies ob served in the infrared and the submillimeter (IDraine et al.l 
2007). In this model, a mixture of dust grains is heated by a 
distribution of starlight intensities. A power law distribution of 
interstellar radiation field intensities, between U m - m and U max , 
represents photo-dissociation regions (PDR), while a constant 
interstellar radiation field {/„„•„ illuminates the diffuse medium 
of the galaxy. 

Different dust mixtures and dust size distributions exist in 
this model, among which three ha ve been tailored to the LMC 
environment and o ne for the SMC (Weingartner & Draine l200ll 
Draine 2002). However, the SMC dust model was tailored 
to SMC extinction curves built on observations of only a few 
stars that mi ght not represent the inter stellar medium of the SMC 
as a whole dGordon & Clavtonlll998l) . Indeed, fitting the "typi- 
cal" SMC extinction curve ( no 2175Abump, steep FUV rise), 
(I Weingartner & Dra ine 2001) implied a lack of very small car- 
bonaceous grains and a predominance of silicate dust grains 
in the SMC, which contrasts with the 8/vm emission, the clear 
PAH features and the lack of silicate features (in absorption 
or in emission) as it is observed in man y regions of the SMC 
(ISandstrom et all 1201 Ot iBot et al.l 12010b . Furthermore, study- 
ing 17 nearby g alaxies observed fr om the near-infrared to the 
sub-millimeter, (IDraine et al.ll2007l) showed that in no case the 
SMC bar dust model was preferred to the Milky-Way (MW) 
or LMC dust model. Given these consider ations, we chose to 
model the dust emission in the SMC with the lDraine & Lil(l2007l) 
dust model using MW type dust. Note however that this choice 
should not significantly affect the far-infrared part of the spec- 
trum. 

The LMC and SMC integrated fluxes from 1.27 to 240/um 
were thus fitted with LMC dust models and MW dust models 
respectively. To do the fit, we use the model spectra as observed 
with the DIRBE, IRAS and Spitzer photometric bands (i.e. tak- 
ing into account the spectral response of the instruments and 
color corrections with respect to the reference spectral shapes) 
and compare them to the observed fluxes in these same bands, 
which avoids the use of color-corrections. The best model dust 
spectra obtained from the fitting procedure are presented in 
Figure [1] and the corresponding parameters obtained are sum- 
marized in Table Q] 

Free-free emission in the Magellanic Clouds is taken from 
the best estimates determined from the best fit of radio data by 
Ilsrael et al.ld20Tob (i.e. S f 2 ^ H { ree = 136 and 12.5Jy for the LMC 
and SMC respectively, with a v~~° A law, and synchrotron spectral 
indexes are qlmc - -0.70 and asMC - -1-09.) The correspond- 
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Fig. 1. Spectral energy distributions from the IR to the radio for the integrated Magellanic clouds. The iDraine & Lil J2007) dust 
model spectra as observed in the different photometric bands (blue squares) is fitted to the observed data points (black stars) from 
1 .27 to 240^m. Synchrotron and free-free spectra are adjusted on the WMAP and radio data. The sum of the three component fit is 
displayed as a dashed line. Error bars correspond to Icr uncertainties. 



Table 1. Results from dust emission modeling of the inte- 
grated spectra: minimum and maximum radiation fields (£/„„■„ 
and U max , in units of the local interstellar radiation field), PDR 
fraction (y), PAH mass fraction {qpAn) and total dust mass (in 
M ) 



Parameter LMC 



SMC SMC (with U m 



0.1) 



model 


LMC 


MW 


MW 




2.5 


2.0 


0.1 


^ max 


10 3 


10 3 


10 3 


y 


0.5% 


7% 


74% 


Ipah 


2.4 


0.47 


0.47 


M dus , 


3.6 • 10 6 


2.9 • 10 5 


1.1 • 10 6 



ing free-free and synchrotron emission for the LMC and SMC 
are shown in Fig.Q] 

The sum of the three fitted components (dust + free-free 
+ synchrotron) is shown as a dashed line in Figure Q] For the 
LMC, the observed emission is well fitted from infrared to mil- 
limeter, as well as in the radio domain. To pHat fluxes at (sub- 
)millimeter wavelength (475/im to 1.22 mm Agui rreet al.ll2003l) 
are well reproduced by the expected dust emission fitted to the 
infrared emission. We observe a small, but significant excess (6- 
7cr) emission in the first three WMAP bands at 3.3, 5 and 7.5 mm 
(90,60 and 40 GHz), above extrapolated models (dust, free-free 
and synchrotron). For the SMC integrated spectrum, a clear ex- 
cess is observed in all TopHat and first four WMAP bands (i.e. 
from 476/im to 10 mm). This shows that in the SMC, the sub- 
mm part of the SED can not be well accounted for with dust 
modelling from the knowledge of the far-infrared peak emission 
(up to 240/im) only. 

Part of the excess observed in the SMC could be due to dust 
colder than can be extrapolated from far-infrared data points. In 
particul ar, a ~ 60/im excess has been reporte d in the Magellanic 
Clouds dBot et al.ll20"oll:lBernard et alJl2008h . This 60/rni excess 
could be a significant part of the total emission in the SMC and 
bias dust models toward high temperatures. To check this, we 
performed a new fit of the SMC spectrum, similar to the mod- 
elling described above, but with a fixed minimum radiation field 
intensity: U m i„ = 0.1 (which corresponds to Tdust ~ and is 




Fig. 2. Same as Fig. Q] right panel, but the minimum radiation 
field in the dust model is fixed to a value of 0. 1 times the solar 
neighborhood value (Tdust ~ 12K). Error bars correspond to lcr 
uncertainties. 



the lowe st radiation field ava ilable in the pre-computed spectra 
given for lDraine & Lil d2007f) dust model). The result of the fit is 
shown in Figure [2] This model better reproduces observed dust 
emission, including TopHat data between 476 /vm and 1.2 mm 
and will be chosen as our best model of dust emission in the 
SMC in order to be conservative. Still, a significant excess emis- 
sion remains above the dust, free-free and synchrotron emission, 
between wavelengths of 3.3 and 10 mm. 

The spectrum of this excess in both galaxies is shown in 
Figure[3]and Table|2] It is obtained by subtracting the best model 
emissions (dust, free-free and synchrotron) to the observed spec- 
tra. The spectral shape and surface brightness of the excess ob- 
served in the LMC and in the SMC are similar (c.f. Figure [3). 
This suggests that they probably arise from the same physical 
process. 
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Table 2. Surface brightness of the observed excess emission in the Magellanic Clouds. Error bars represent 1 <x uncertainties. 
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Q. 




MJy/sr 


MJy/sr 


MJy/sr 


MJy/sr 


MJy/sr 


MJy/sr 


sr 


SMC 


0.019 ±0.014 


0.018 ± 0.002 


0.0075 + 0.0011 


0.0035 ± 0.0006 


0.0024 ± 0.0004 


0.0014 ± 0.0002 


0.00544 


LMC 


< 0.032 


0.012 ± 0.001 


0.0040 ± 0.0005 


0.0014 ± 0.0002 


0.0007 ± 0.0002 


< 0.0003 


0.0160 



150.0 
0.100 



75.0 



v (GHz) 



50.0 



37.5 



0.010 



0.001 



SMC 

LMC 



dB(2.725K)/dT.AT 
cold dust fit for the LMC 
cold dust fit for the SMC 



2.0x10 3 4.0x10 3 6.0x10 3 8.0x10 3 1.0x10" 1.2x10" 1.4x10* 

Fig. 3. Spectral energy distribution of the mm-cm excess in the 
LMC (in green) and in the SMC (in black). The spectral shape 
of th e CMB fluctuations (d B y (2. ,7 '25 K)/dT x DeltaT iFixsenl 
2009) is scaled to the SMC excess brightnesses for comparison. 
Modified black body fits to the LMC and the SMC excess rep- 
resenting a possible cold dust component are shown as dotted 
lines. 



3. CMB fluctuations? 

The intensity of the cosmic microwave background fluctuations 
are expected to be non-negligible at low surface brightnesses 
like the one observed for the Magellanic Clouds in the mm-cm 
range. Moreover, in the SMC, the spectral energy distribution 
of the excess is of similar shape to the CMB fluctuations (c.f. 
Figure[3]). We therefore estimate the probability that CMB fluc- 
tuations in the background of the Magellanic Clouds are at the 
origin of this submm-cm excess. 

To do so, we created a CMB simulati on map that repro - 
duces the observed CMB power spectrum dLarson et al. 2010). 
We draw 5000 random positions in this CMB map and extract 
intensities of CMB fluctuations at 5mm (60GHz) inside an area 
equal to the one of the SMC and with an annular background 
subtraction. In Figure |H we plot an histogram of CMB fluctu- 
ations intensities obtained for these 5000 random regions. We 
then estimate the probability of observing CMB fluctuations at 
a surface brightness that is compatible with the observed excess 
surface brightness (c.f. Tab. [2]) at 5mm within 3<x . We find that 
this probability amounts to 5% for SMC. 

We did the same estimate for the LMC. The histogram of 
CMB fluctuations intensities obtained with the same method is 
shown in Figure [2] (left panel). The difference observed between 
the histograms for the LMC and the SMC corresponds to a size 
difference between the regions considered in each case (~ 4° 
and ~ 2° radii respectively). For the LMC, we find a 14% prob- 
ability that the intensity of the excess at 5mm is due to CMB 
fluctuations. However, the spectral shape of the excess in this 
case differs from the one of CMB fluctuations. 



Furthermore, the fact that both Magellanic Clouds show a 
similar mm-cm excess emission despite their different size on 
the sky argues against the hypothesis that this excess is due to 
CMB fluctuations. 

We therefore consider it very unlikely that the excess in 
the Magellanic Clouds originates from CMB fluctuations in the 
background of these galaxies. We note that our test is based on 
average surface brightnesses and as such does not make any use 
of the spatial distribution of the excess in the Magellanic Clouds. 
A study of the excess with respect to the expected CMB fluctua- 
tions at different spatial scales would help to better discriminate 
the significance of CMB fluctuations in that respect. 

4. Very cold dust? 

It is tantalizing to try to explain the excess observed with very 
cold dust, especially since colder dust in the SMC can ex- 
plain the sub-millimeter part of the excess (c.f. Section [2j. To 
check this hypothesis, the excess observed is fitted by a mod- 
ified blackbody representing a possible cold dust component. 
Best fits are shown in Fig.[3]and give dust temperatures of ~ 3K 
and spectral indices of B<j°$ c = 1.2 and Bf$ c = 0.33. Such a 
cold thermal equilibrium temperature for big grains is unreal- 
istic. Furthermore, the spectral dust emissivity index deduced 
for the SMC is also unrealistically low. Indeed, current obser- 
vations suggest R to be in the 



e range l-2.5(Boulan geret al.fi 996: 
et al.l2003aHParadis et al . 2009) and 



Laga che et alJ 19991: iDupac i 
the Krame rs-Konig equat ions suggest that a value of 1 as a lower 
limit to B (Emerson 1988). The existence of very cold dust that 
would explain the observed millimeter excess therefore seems 
very unlikely. 

5. A change of the dust emissivity spectral index? 

Another hypothesis for explaining the millimeter excess is a 
change of the spectral emissivity index from the F IR to the mm. 
Such a change has been observed in our Galaxy (Paradi s et alJ 
120091) and is related to intrinsic dust pro perties. To test thi s 
assumption, we used a model proposed by iMenv et al.l (|2007), 
based on the properties of amorphous dust grains. They consider 
two complementary physical models. The first one describes 
excitation of acoustic lattice vibrations, due to the interaction 
between electromagnetic fields and a disordered charge distri- 
bution (DCD). This DCD model is characterized by a correla- 
tion length, used to describe the disorder degree of amorphous 
state. The second model focuses on the processes associated 
with a distribution of low-energy two-level systems (TLS), such 
as resonant absorption and relaxation processes, temperature- 
dependent emission. This DCD/TLS model was succesfully used 
to reprod uce the FIRAS/WMAP spectrum of gala ctic diffuse 
emission dParadisll2007blParadis et al1l2010. in prej2j). 

Fits to the LMC and SMC integrated spectra were done for 
all points from 100^m to 9.1mm. To do so, three parameters 
were set to be free: the dust temperature, the correlation length 
(L c ) and the ratio of TLS over DCD effects (A ratio). We did 
not attempt to fit the 60//m data point as the emission from very 
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Fig. 4. Histogram of the simulated intensity fluctuations of the CMB measured similarly to the fluxes of the LMC and SMC (left and 
right panels, respectively) at 5mm (60 GHz), for 5000 random positions on the sky. The histograms have been normalized to their 
peak value and the bin size is 5 x 10~ 4 MJy/sr). The solid vertical line in each pannel corresponds to the measured excess intensity 
at this wavelength, while the dashed lines give the 3<x error estimates to either side of this measured excess. 



small grains that are transiently heated is known to be signif- 
icant at that wavelength. Best fits are shown in Figure [5] and 
correspond to dust equilibrium temperatures of 2 IK and 24K, 
correlation lengths of 6.4 and 12.85 nm and A ratios of 4.4 and 
19.4, for the LMC and the SMC respectively. 

In the LMC, the DCD/TLS model is able to reproduce ob- 
servations, assuming a decrease of the correlation length by a 
factor of 2 and an increase of the TLS processes intensity by a 
factor of 1.8, with respect to the parameters used to reproduce 
the spectrum of galactic diffuse emission. This result indicates a 
good agreement between dust properties in the diffuse medium 
of our Galaxy and in the overall LMC, even if the amorphization 
degre43 of dust grains is different. 

However, for the SMC, the flattening of the spectrum is too 
pronounced to be explained only by intrinsic dust properties. 
Figure|5]shows a fit of the SMC SED with the DCD/TLS model. 
In this case, TLS effect intensity had to be enhanced by a factor 
of 8. Despite this enhancement, an additional component in the 
millimeter domain is needed to reproduce the data. 

Hence, taking into account the properties of amorphous 
grains satisfactorily explains the excess in the LMC but the 
model fails at reproducing the one observed in the SMC. 

6. Spinning dust emission? 

The main mechanism invoked to explain Galactic anomalous ex- 
cess at milli meter wavelengths is spinning dust emission, as first 
proposed by Draine & Lazarianl d!998bl) . Unlike the excess in 
the Magellanic clouds, Galactic excess peaks between 20 and 40 
GHz. However, models predict that the peak frequency of spin- 
ning dust emission shifts with radiation field intensity, interstel- 
lar gas physical parameter fl grain size dis t ribution and electric 
dipol e moment dAli-Haimoud et al] 2009; lYsard & Verstraetd 
l2010h . 

Here, we explore the possibility that the excess observed 
in the Magellanic Clouds originates from spinning dust emis- 

1 The amorphization degree indicates if the internal structure of a 
grain is fully or partially amorphous. 

2 Collision with interstellar g as species may excite or damp the rota- 
tion iDraine & Lazarianlll998bh . 



Table 3. Gas parameters obtained with CLOUDY. Densities are 
expressed in cm -3 . 



medium 


Go 


n« 


T 

S as 


n c + 




n H + 




LMC diffuse 


2.5 


30 


189.7 


2.0 x 10 


-3 


6.8 x 10" 


-3 


LMC PDR 


2.5 


10 4 


109.9 


0.24 




0.75 






1000 


10 4 


288.8 


0.66 




1.8 




SMC diffuse 


0.1 


30 


124.1 


6.5 x 10" 


-4 


3.6 x 10" 


3 


SMC 


0.1 


10 J 


115.2 


3.7 x 10" 


-3 


3.8 x 10" 


-2 


moderate PDR 


1 


10 3 


127.7 


1.9 x 10" 


_2 


8.0 x 10" 


-2 




1000 


10 3 


1168.2 


2.2 x 10" 


-2 


0.88 




SMC 


0.1 


10 5 


63.2 


3.6 x 10" 


-4 


0.3 




dense PDR 


1 


10 5 


124.0 


2.1 x 10" 


-2 


1.4 






1000 


10 5 


278.7 


2.2 




9.3 





sion produced by very small dust grains (i.e. PAHs). To do so, 
we fit the excess observed in the Magellanic Clouds with a 
model of spinning dust. In both cases, and in order to reduce 
the number of free parameters, we make some ad-hoc assump- 
tions. First, the electric dip ole moment distribution is yU = m yN 
dDraine & Laza rian 1998b), where N is the number of atoms in 
the gr ain and m is a constant equal to 0.4 D dYsard & Verstraetel 
2010). The size distribution of dust grains is taken to be the 
same as in the iDraine & Lil (12007) model (in order to be consis- 
tent with the fit of the mid-infrared part of the SEDs). We make 
the assumption that the interstellar medium of the Magellanic 
Clouds can be modelled as the sum of a diffuse medium and a 
PDR phase. The diffuse medium is assumed to have the same gas 
density as in the cold neutral medium in our Galaxy (n# ~ 30 
cm -3 ), while the PDR component has the same density as the 
Orion Bar (n# ~ 10 4 cm" 3 ). Using the radiation field distri- 
bution described in Tab. [1] other relevant gas p arameters are 
calcu lated at thermal equilibrium with CLOUDY (Ferlan d~et al.l 
1998). We assume Z$mc - Z Q /6 and Zimc — Z Q /2. Parameters 
were taken from the optically thin zone of isochoric simulations 
with CLOUDY. The main gas species abundances vary with ra- 
diation field and density. We present some of the gas parameters 
obtained in Table [3] for the different phases in each galaxy and 
for a sample of radiation field intensities. 
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Fig. 5. SEDs of the Magellanic Clouds in the FIR to mm regime, once the free-free and synchrotron emission removed. This spectral 
energy distribution is fitted from lOOyum to 9.1mm with a dust model with DCD/TLS effects. The spectral shape of the DCD effect 
is shown in light blue while the TLS effects are displayed in dark blue and the total emission is displayed in red 



For the LMC, a good fit is achieved and is shown in Fig. [6] 
(left panel) overlaid on the observed excess emission. The re- 
sulting PAH mass fraction is qpAH =2.15%, with a carbon abun- 
dance in PAH of 0.9 ppm in the diffuse medium and 3.06 ppm in 
the PDR. 

For the SMC however, no good fit was achieved with a dif- 
fuse+PDR model. Separating the PDR component into a dense 
PDR component (with n# = 10 cm -3 ) and a moderately dense 
PDR component (with Hh = 10 3 cm~ 3 ), a good fit of the excess 
emission observed in the SMC is achieved and is shown in Fig. 
|6] (right panel). In this case, we deduce qpAH = 0.39%, with a 
carbon abundance in PAH of 2.20 ppm in the diffuse medium 
and 0.49 ppm in the dense PDR and 0.47 ppm in the moderately 
dense PDR. 

The spinning dust emission model can therefore reproduce 
the excess observed in both Magellanic clouds with parameters 
that are consistent with our current understanding of the inter- 
stellar medium in these environments. In particular, the PAH 
mass fractions obtained from the best fits of the millimeter ex- 
cess are consistent with those obtained from the far-infrared dust 
emission modelling of the spectral energy distributions. This is 
consistent with the idea that spinning dust emission is produced 
by the smallest grains which are responsible for the mid-IR emis- 
sion. 

If the millimeter-centimeter excess is due to spinning dust, 
the peak frequency of spinning dust emission in the Magellanic 
Clouds (139 and 160 GHz - 2.2mm and 1.9mm- for the LMC 
and SMC respectively) is shifted with respect to what is ob- 
served in our Galaxy. This shift depends on three main parame- 
ters: the size (the smaller the grain, the higher its emission fre- 
quency), the density (if the medium is denser, more grain-gas 
collisions occur and the anomalous emission peak is shifted to- 
ward higher frequencies) and the intensity of the Interstellar ra- 
diation Field (ISRF) (above a threshold of G = 10, if the ISRF 
intensity increases, the emission peak is shifted to higher fre- 
quencies). In order to reproduce the millimeter excess observed 
in the Magellanic Clouds, all three parameters were equally im- 
portant. 

This study shows that spinning dust emission is a possi- 
ble solution for the nature of the observed excess, in the two 
Magellanic Clouds. However, without complementary informa- 



tions, many assumptions had to be made (e.g. mixing of two 
or three components along the line of si ght, density of the 
interstellar medium ). Also, recent studies dHoang et al.l 12010: 
Silsb ee et al.l 12010) have shown that refinements on spinning 
dust models can lead to significant changes in the shape and in 
intensity of the spinning dust emission peak. The inclusion of 
such effects are beyond the scope of this paper, but could help the 
understanding of the excess in the SMC in terms of spinning dust 
emission. We do not claim that the best fits obtained are a unique 
solution. Higher resolution studies with Planck and Herschel 
will be necessary to confirm this origin of the excess by assess- 
ing the distribution of the millimeter excess in the Magellanic 
Clouds as a function of local conditions in better defined envi- 
ronments (e.g. diffuse medium only, HII regions). For example, 
it would be interesting to try to correlate the mm/cm excess in 
the Magellanic Clouds with PAHs emission to test the scenario 
of spinning dust emission. 



7. Conclusion 

Comprehensive spectral energy distributions of the Magellanic 
Clouds from the near-infrared to the radio are modelled with 
dust, free-free and synchrotron emission. An excess above these 
three components is observed at millimeter wavelengths. This 
excess has a similar spectral shape in both Magellanic Clouds 
but is more prominent in the SMC total SED than in the LMC. 

We explore different scenarios for the origin of this excess 
emission in the millimeter regime: very cold dust, CMB fluctu- 
ations, a change of the dust emissivity spectral index and spin- 
ning dust emission. It is shown that very cold dust grains imply 
a thouroughly unrealistic nature for this excess and CMB fluctu- 
ations are unlikely to solely create such an excess. 

For the LMC, the mm-cm excess can be explained equally 
well with TLS/DCD effects in amorphous grains or with spin- 
ning dust emission. 

This is drastically different in the SMC, where the mm-cm 
excess is more pronounced. In this case, CMB fluctuations be- 
come very unlikely and the inclusion of TLS/DCD effects alone 
does not reproduce the excess. However, spinning dust emission 
models are capable to explain the observed excess with parame- 
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Fig. 6. Best fits of the millimeter excess observed in the Magellanic Clouds with anomalous dust components as described in the 
text. 



ters that are consistent with our understanding of the ISM in the 
SMC, but many assumptions had to be made to obtain this result. 

We hence confirm the existence of an unexpected mm-cm 
excess in the Magellanic Clouds, but the nature of this excess 
remains unclear. Spinning dust is a promising solution, but this 
should be further tested in better defined physical environments. 
The mm-cm excess observed in the Magellanic Clouds could 
easily be of a more complex nature, having multiple causes 
rather than a single one, so that the different effects explored 
here might add up. In this context, Planck and Herschel observa- 
tions of the Magellanic Clouds are needed to map this excess at 
a better spatial resolution, trace its spatial variations and probe 
its origin more definitely. 
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